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Abstract 

We study the artificial molecular states formed in laterally coupled double semiconductor nanor- 
ings by systems containing one, two and three electrons. An interplay of the interring tunneling 
and the electron-electron interaction is described and its consequences for the magnetization and 
charging properties of the system are determined. It is shown that both the magnetic dipole mo- 
ment generated by the double ring structure and the chemical potential of the system as function 
of the external magnetic field strongly depend on the number of electrons and the interring barrier 
thickness. Both the magnetization and chemical potentials exhibit cusps at the magnetic fields 
inducing ground-state parity and / or spin transformations. The symmetry transformations are 
discussed for various tunnel coupling strengths: from rings coupled only electrostatically to the 
limit of coalesced rings. We find that in the ground-states for rings of different radii the magnetic 
field transfers the electron charge from one ring to the other. The calculations are performed with 
the configuration interaction method based on an approach of Gaussian functions centered on a 
rectangular array of points covering the studied structure. Electron-electron correlation is also 
discussed. 

PACS numbers: 73.40. Gk 
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I. INTRODUCTION 



In the eighties semiconductor rings of micrometer size were investigated in quest for sig- 
natures of the Aharonov-Bohm effect^ in conductance measurements.-*^ In the next decade 
magnetization produced by persistent currents circulating around semiconductor rings was 
measured.- Subsequent technological advances allowed for fabrication of rings with nanome- 
ter radii with detectable quantum size effects. Nowadays, the quantum rings are produced 
with the etching^ or surface oxidation techniques^ as well as grown by self-assembly.^ Trans- 
port experiments are performed on open quantum rings,- while the closed rings are studied in 
the context of single-electron charging'' or optical properties.->i2iiii Recently, magnetization 
signal of large ensembles of semiconductor self-assembled nanorings has been detected.— 

The theoretical literature on quantum rings is very rich. The authors mainly con- 
centrated on properties of a single isolated quantum ring.i2iiiii2iiiiii^ii^iiZii^ii£ At present 
there is a growing interest in systems of multiple quantum rings including arrays of quan- 
tum rings.— 1^ Moreover, double rings are produced in both concentrio^^i^^i^i and vertical^ 
configurations. The tunnel and electrostatic coupling was theoretically studied for both 
concentrio^^iSIiS^iSSiSSiSiiS^i^^ and vertically stacked rings.— Recently, states of a single-electron 
in a pair of laterally coupled quantum rings were described.— The magnetization generated 
by planar arrays of interacting quantum rings^^ for neglected tunneling between the rings 
was also discussed. The purpose of the present paper is to describe the system of up to three 
electrons in an artificial molecule formed by two quantum rings with the account taken for 
both the tunnel coupling and the electron-electron interaction. We investigate the compe- 
tition between the tunnel and Coulomb coupling, the electron-electron correlation as well 
as the charging and magnetization properties. The numerical results are provided for the 
etched InGaAs/GaAs rings^ with low indium concentration and consequently low potential 
depth which favors electron tunneling between the rings. 

The evolution of the single-electron ground-state with the magnetic field for artificial 
molecules formed by quantum rings is significantly more complex than for the double quan- 
tum dots. In double dots the role of the magnetic field for the eigenstates of the single 
electron is limited to reduction of the interdot tunnel coupling.— In double rings the mag- 
netic field drives the angular momentum transitions of the single-electron within each of 
the rings.-i^ When the tunnel coupling is activated the single-electron ground-state becomes 
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localized at the contact between the rings^ with the strength of the localization oscillating 
in function of the magnetic field with a period corresponding to the flux quantum thread- 
ing the rings. For identical rings the single-electron ground state corresponds to a binding 
orbitaF^ and possesses an even spatial parity irrespective of the value of the magnetic field. 
In this paper we show that both the spin and spatial symmetry transitions occur with the 
external magnetic field^^ double rings containing few electrons. We find that the symme- 
try transformations in the external magnetic field depend strongly on the interring barrier 
thickness and they have quite a different character for various numbers of electrons confined 
within the double ring structure. Since the transformations influence strongly the charging 
and magnetic properties, an evidence of the tunnel coupling between the rings should be 
detectable by measurements of chemical potentials and magnetization. 

Our discussion covers asymmetric configurations composed of rings of different radii for 
which we find that the magnetic field induces an oscillatory switching of the ground-state 
localization from one ring to the other. This effect can be used to transfer the electron 
between the rings without a need for applying an external electric field to the system. 

In this work we use the configuration interaction approach which allows for a numerically 
exact solution of the Schroedinger equation for a few electrons. Application of the config- 
uration interaction method to the laterally coupled double ring structure in the external 
magnetic structure is challenging as compared to both laterally coupled dots and double 
rings in concentric and vertical configurations. The basis used for the configuration inter- 
action calculation has to keep track of the angular momentum transitions which occur in 
each of the rings separately but the total angular momentum cannot be used as a quantum 
number for selection of the basis set since the system does not possess circular symmetry 
unlike concentric double or vertically stacked rings. For the purpose of the present study we 
developed quite a powerful technique in which a Gaussian functions are used with centers 
distributed on a regular array. The presented technique is universal and can be applied to 
few-electron systems in arbitrary smooth confinement potentials. 

This paper is organized as follows. In Section II we present the model Hamiltonian, 
the confinement potential and the configuration interaction approach based on a mesh of 
Gaussian functions. A test of the method for two-electron states is presented and the 
limitations of the approach are explained on this example. The results are given in Section 
III. We start by single-electron states of a single and identical double rings. The discussion is 
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then extended to few-electron states and non-identical rings. The summary and conclusions 
are given in Section IV. 



II. THEORY 

We consider the following few-electron Hamiltonian 

N N 1 

H = j:h^+ E — (1) 

i=i j=i,j>j^'«i 

where hi is single electron energy operator and apply the configuration interaction approach 
in which H operator is diagonalized in a basis of many electron wavefunctions with deter- 
mined values of total spin S and its projection on z axis 5*^. The basis functions are generated 
with the help of projection operator^'^ as the linear combinations of Slater determinants built 
of eigenfunctions of a single-electron Hamiltonian: 

We use the vector potential in a symmetric gauge: A(r) = B/2{—y, x, 0). The single-electron 
eigenproblem is diagonalized in a basis 

N 



m = Y.c^^f^^ (3) 

where the basis function have the form 

fair) = exp ( - exp ( - ^(B X R„) ■ r) . (4) 

In Eq. m the probability density associated with each basis function is a Gaussian centered 
at point Rq = {xa,ya), cr is responsible for the strength of the localization and the term 
with the imaginary exponent introduces the magnetic translation, which ensured the gauge 
invariance of the basis, or in other word equivalence of all the basis functions irrespective of 
the localization center. The centers are distributed on a regular mesh of points (see below). 
The matrix elements of the Coulomb interaction are integrated according to the procedure 
explained in the Appendix. 

We model the potential of a single two-dimensional quantum ring by the formula: 



where Ri(r) denotes the center of the left (right) ring. We assume a = 20 for which the 
potential is nearly a square quantum well. The other parameters of the rings are adopted 
for the etched Ino.1Gao.9As/GaAs quantum rings^i^^: Vq = 50 meV, the radius of single 
quantum ring Rq = 30 nm, (Tq = 20 nm, effective mass of an electron m* = 0.05 and 
dielectric constant e = 12.4. 

Confinement potential of the double ring structure is assumed in form, 

Ve.t{r)=mm(Vi{r),Vr{r)^ (6) 

where Vi{r),Vr{r) are the confinement potentials of the left and the right quantum ring 
respectively. We use a dimensionless parameter d which describes the distance between the 
centers of two coupled rings: 

d=\^i^ (7) 



Ro 

The confinement potentials of two laterally coupled quantum rings for various d are shown 
in Fig. [H 

The centers of the single electron basis functions (jl]) are distributed on a mesh of x ky 
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FIG. 1: Double ring structures for various parameters of the interring distance d and the confine- 
ment potential cross sections through the centers of the rings (b) for three values of parameter 
d. 

points, where k^ = ^^'l'^'^^ x ky, ky = 40. For all functions fa value of a parameter is identical 
and equals a = yff^- This value was optimized for the description of the low-energy part 
of the single-electron spectrum. 
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We have tested our method for the problem of two electrons confined in two-dimensional 
harmonic oscillator potential, which can be easily solved with the center-of-mass separation 
technique with an arbitrary precision (the results can be treated as "exact"). We assumed 
the oscillator energy huj = 1 meV and the calculations were performed in the square region of 
size 20/0 X 20/0 where /q = \jTi/m*uj is the oscillator length. We have solved the eigenproblem 
with the configuration interaction method for basis containing = 20 x 20 and = 50 x 50 
centers. The exact energy spectra and the results of the present approach are plotted in 
function of B in Fig. [2j 




FIG. 2: Comparison of the exact two-electron energy spectrum (solid lines) for the 2D harmonic 
oscillator confinement potential and the results obtained in the present CI method (red dots) for 
the mesh of 50 x 50 and 20 x 20 Gaussians. The energy spectra are calculated with respect to 
the non-interacting ground state energy equal to double of the Fock-Darwin^E ground-state energy 
level. 

For the basis with = 50 x 50 elements, the lower part of the energy spectrum agrees 
very well with the exact results for magnetic fields up to 4 T. However, for basis with 

= 20 X 20 elements, the energy of the ground state is visibly overestimated and the 
overestimate becomes quite large above 3 T. 

For both basis the probability densities are circularly symmetric about the center of the 
dot for 5 = 1 T and = 3 T (see Fig. [3](a,b) and [3];d,e)). For 5 = 4 T only the 
probability density calculated for A^ = 50 x 50 elements still preserves the confinement 
potential symmetry. In the right panel of Fig. [3] we plot the mean value of the angular 
momentum for the smaller basis. With the magnetic field the ground-state becomes more 
strongly localized and acquires a high angular momentum. For fixed number of centers there 
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FIG. 3: (a-f) The total probability densities for = 50 x 50 and = 20 x 20 basis centers 
for two interacting electrons confined in a two-dimensional harmonic oscillator potential {huj = 1 
meV) in the singlet states. The dots on (c) and (f) indicate the centers of Gaussian functions, (g) 
Expectation values the total angular momentum for 20 x 20 elements (points) and 50 x 50 Gaussian 
(crosses). 

is a limitation to the maximal angular momentum that can be described in the single-electron 
basis. The single-electron states of high angular momenta correspond to sign oscillations of 
both the real and imaginary part over the angular variable. Clearly, the number of centers 
in the region where the probability density is non-zero sets the limitation to the maximal 
frequency of the oscillations. For B = 4 T the angular momentum obtained in the basis 
of 20 X 20 elements is non-integer (L^) = 6.86 and smaller than the exact ground-state 
value of (Lz) = 8, which is nevertheless reproduced by the 50 x 50 basis. Concluding, the 
present approach allows for a nearly exact solution of the few-electron Schroedinger equation. 
There is a limitation to the maximal angular momentum that can be accounted for, but the 
accuracy of the results is easily verified by comparing results of the method produced by 
various meshes. 



III. RESULTS 

A. Single-electron single-ring states 

The single-electron spectrum for a single quantum ring with potential ([6]) is shown in Fig. 
m The energy levels plotted with black lines correspond to the lowest state of the radial 
quantization and the red ones to the first radial excitation. The energy spacing between 
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these two branches is only about 10 meV, which results of the small depth of the studied 
structure. 

Numbers in Figure H] denote the angular momentum. In the ground-state the magnetic 
"period" of the ground-state angular momentum transition is -Bp = 1.05 T which corresponds 
to the flux quantum [BpirR^ = e/h) for a strictly one-dimensional ring of radius -Ri^ = 35.4 
nm, slightly larger than the Rq parameter (assumed equal to 30 nm), and closer to the 
average distance from the center of the ring 37 nm calculated over the radial coordinate 
of the ring. The magnetic period for the first radial excitation is much larger and equals 
Bp = 1.5 T, which corresponds to the ID ring of Rid = 28.9 nm. The smaller effective 
radius for the excited radial state than for the ground-state may be surprising since usually 
the excited states occupy larger area than the ground-state. In order to explain this feature 
we plotted the electron densities in Fig. [5] for zero magnetic field. Fig. [5](a) shows the 
density for the ground-state and Fig. Mjo) the density for the L = state of the first radial 
excitation. Fig. [5] (c) and (d) show the densities of the angular momentum L = 3 eigenstates 
- the lowest-energy state and the first excited state, respectively. We notice that the radial 
wave functions for L = are more strongly localized than the ones for L = 3. This is due 
to the centrifugal effective potential present in the single-electron Hamiltonian for a 
single ring written in cylindrical coordinates for the angular momentum L eigenstate 

R'S' '''i' ,„ > 2 1. r 

where p = y/x^^flp- Fig. [5] shows that states that correspond to the radial excitation 
with electron density forming two concentric rings occupy indeed a larger area than those 
corresponding to the lowest radial state. However, the inner density ring contains most of 
the electron charge and its radius is smaller then the electron density ring in the lowest- 
energy radial state [cf. Fig. [S](a) and (b), as well as Fig. MyC) and (d)], which explains why 
the effective radius value Rid obtained for the branch of the excited energy levels is smaller 
than for the ground state. 

B. Single electron states for a pair of identical rings 

The coupling between the rings is of a pure tunnel character for the single-electron states. 
The single-electron energy spectra are given in Fig. [6] for various interring distances. For 
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FIG. 4: Single-electron spectra for a single quantum ring. Black curves show the energy levels 
associated with the lowest radial state, and the red curves the energy levels related to the first 
radial excitation. Numbers in the figure indicate the angular momentum. 
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FIG. 5: Electron density for the single electron in a single ring in zero magnetic field. Plots (a,b) 
correspond to zero angular momentum and (c,d) to L = 3. Densities in (a,c) are obtained for the 
lowest-energy radial state, and (c,d) for the first radial excitation. 

d = 3.75 the tunnel coupling is negligible and the spectrum [see Fig. [6](a)] consists of energy 
levels which are twofold degenerate with respect to the parity. This energy spectrum is 
in fact identical with the one of a single separate ring (see Fig. H]). In Fig. [6](b) plotted 
for d = 3.5, we spot lifting of the even-odd degeneracy which occurs when the interring 
tunnel coupling is activated. Note that when the degeneracy with respect to the parity 
is lifted, lower energy level always corresponds to the even parity energy level, i.e. to the 
binding orbital. The even-odd degeneracy lifted by the tunnel coupling is restored for higher 
magnetic fields. This results of the attenuation of the tunnel coupling by the magnetic field 
which enhances electron localization within each of the rings. For d = 3.25 - when the rings 



nearly touch one another - the even and odd energy levels differ significantly in the energy. 
We observe a pronounced avoided crossings in the spectrum that occur separately for the 
even and odd energy levels. The pattern of the avoided crossings in the even and odd parts 
of the spectrum is similar, only the avoided crossings for the odd energy levels are narrower. 
The odd parity energy levels correspond to wave functions vanishing in the center of the 
tunnel barrier so the tunnel coupling between the rings is naturally less pronounced than 
for the even energy levels. 

In the absence of the tunnel coupling between the rings {d > 3.75) the electron density 
of the stationary states reproduces the circular symmetry of separate rings - irrespective of 
the magnetic field value. A distinct dependence on the magnetic field occurs only when the 
tunnel coupling is present. The charge densities for the lowest even and odd energy levels 
are given for d = 3.25 in the lower (a-d) and upper (e-h) rows of Fig. [71 respectively. The 
lowest even (odd) energy level is the ground state (the first excited state). In the ground 
state the electron tends to stay near the symmetry center of the double ring system. The 
extent of its localization varies with the magnetic field. The strongest localization near the 
interring contact area [cf. Fig. [7](c)] corresponds to the center of the avoided crossing of 
the two lowest energy levels [see the spectrum of Fig. Et^c)]. In the odd parity energy level 
the electron is by the symmetry forbidden to be found at the center of the structure and it 
tends to occupy the extreme ends of rings. The strength of localization at the ends oscillates 
with the magnetic field and is the largest near the center of the avoided crossings occurring 
between two lowest odd energy levels [see Fig. [71(g) for B = 1 T and Fig. [6](c)]. 

Coalesced rings at d = 3 form a cavity of increased width near x = (see Fig. 1), 
which acts like a quantum dot - in the ground-state the electron becomes localized at the 
interring contact area [Fig. [8](a-c)] and only weakly penetrates the more distant parts of 
the structure. The ground-state electron density weakly depends on the magnetic field since 
the electron is localized within an area which is quite small. The energy of this (ground) 
state is distinctly lowered with respect to the excited part of the spectrum [compare Fig. 
M^d) for = 3 Fig. [6|(c) and d = 3.25]. The excited energy levels are essentially unchanged 
with respect to d = 3.25, only the avoided crossings become larger. Fig. [8] shows that the 
densities in the first excited state of the coalesced system are similar to the d = 3.25 case: 
they stay spread all over the double ring structure and are not limited to the contact area 
as in the ground-state. In the spectrum calculated for d = 3.25 [Fig. [6](c)] a trace of the 
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Aharonov-Bohm oscillation for the ground-state energy is still visible since at least outside 
the center of the avoided crossings the ground state wave function encircles the rings (see 
Fig. [7]). For d = 3 even this residual oscillation of the ground-state energy disappears [see 
Fig. E^d)]. However, oscillations are still present in the excited part of the spectrum in 
which the wave function covers the entire structure. Therefore the coalescing of the rings 
mainly perturbs the ground state which becomes localized in the quantum dot formed at 
the contact of the rings and not the excited states which are spread all over the double ring 
structure. 
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FIG. 6: Single-electron spectra for a pair of identical rings for various interring distances. Solid 
lines show the even parity energy levels and the dashed lines the odd parity energy levels. In (a) 
the angular momenta with respect to the center of the ring are listed. 




FIG. 7: Lower row of plots (a-d) shows the charge density in the ground-state (of the even parity) 
for d = 3.25 (the energy levels are given in Fig. EJc)). The upper row of plots shows the first 
excited state, which is the lowest-energy state of the odd parity. Columns correspond to different 
magnetic fields. 
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FIG. 8: Same as Fig. [7] but for d = 3. The corresponding spectrum is displayed in Fig. [Hd). 
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FIG. 9: Two-electron spectra for a pair of identical rings. Solid and dashed curves show the even 
and odd parity energy levels. Black and red lines correspond to the total spin S = and 5 = 1. 

C. Identical rings: two electrons 

The electrostatic coupling between the rings appears with the second electron introduced 
to the system. The two-electron spectrum for a pair of rings is displayed in Fig. [91 For 
two electrons the Coulomb repulsion makes the interring barrier higher reducing the tunnel 
coupling and in the low-energy part spectrum the electrons tend to occupy separate rings. 
For wide interring barrier and vanishing tunnel coupling the electron separation is complete. 
In consequence , for ci = 3.75, 3.5 and 3.25 the ground state is two-fold degenerate [see Fig. 
ini^a-c)]. One of the ground states is the spin singlet of the even parity and the other is the 
spin triplet of the odd parity. In the two-electron systems the singlet-triplet degeneracy 
(vanishing exchange interaction's,^ occurs when the electrons occupy distinct and separated 
locations. For double dots this degeneracy is found for wide^^ interring barriers or high 
magnetic fields. Similar is the effect of the Wigner crystallization in elongated, quasi one- 
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FIG. 10: Two-electron densities for d = 3.25 for four lowest-energy singlet states (for triplets the 
results are nearly identical). The lower the plot the lower the energy. The ground state (a-c) and 
the third excited state are of the even parity for all B. For B = 1.03 T the first excited state 
is of the even parity and the second excited state is odd. For the two other magnetic field the 
symmetries of this two energy levels are inverted - see also Fig. [9)[c). 
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FIG. 11: Two-electron densities for d = 3 for four lowest-energy singlet states. The lower the plot 
the lower the energy. Plots (a-c) correspond to the ground-state. 

dimensional quantum dots.— 

For a single-electron the effects of the tunnel coupling - splitting of the odd and even 
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parity energy levels and avoided crossings in the odd and even parts of the spectrum - 
were visible already for d = 3.5 [Fig. [6](b)] and for d = 3.25 [Fig. [6](c)] they were already 
quite strong. Lifting of the spin degeneracy is only visible for the coalesced rings [d = 3 
see Fig[n](d)]. Note that also in this case the singlet and triplet becomes degenerate when 
the magnetic field is switched on enhancing the electron localization in opposite ends of the 
structure. 

For wide interring barrier, when the tunnel coupling is absent, the separate rings are 
only coupled by the Coulomb interaction of confined electrons. The electrostatic potential 
of the electron localized in one of the rings perturbs the symmetry of the other ring. In 
consequence the angular momentum of a single ring is not a good quantum number and one 
obtains avoided crossing between energy levels [Fig. [n](a-c)] instead of the crossings [compare 
Fig. [6](a)]. Note, that the center of the avoided crossing for two electrons [see Fig. 2 [9](a)] 
occurs for the value of the magnetic field for which the ground state angular momentum 
transition occurs for a single electron [cf. Fig. [6](a)]. 

Electron densities are displayed in Fig. [10] for d = 3.25. In the ground-state the electrons 
stay at the opposite ends of the rings. The most pronounced separation occurs at the center 
of the avoided crossing near B = 1.03 T. Oscillation as pronounced as in the ground state is 
observed for the third excited singlet state. In this state the electrons are mostly localized 
in the center of the structure in contrary to the ground-state, for which electron mostly 
occupy ends of the double structure. The electron density in the third excited state is in 
fact similar to the single-electron ground state [cf. Fig. [T0](j-1) and Fig. El^b-d)]. However, 
the Coulomb repulsion makes the two-electron density less strongly localized. In the first 
two excited singlets [Fig. [TOT d-i)] the reaction of the density to the field is much weaker 
than for the ground-state and the third excited singlet. 

The singlet electron densities for coalesced rings {d = 3) are displayed in Fig. [TTl Com- 
parison of this plot with the single-electron densities of Fig. [8] indicates that in the two 
lowest-energy two-electron levels [Fig. [TT](a-f)] the single-electron ground-state [with den- 
sity localized in the interring contact area - Fig. fTTT a-c)] and the first excited state [electron 
density localized at the extreme ends of both the rings - Fig. [TTT d-f)] contribute nearly 
equally. In the second [third] excited state the single-electron first excited [ground] state 
contribution is dominant - see Fig. [TlT g-i) [Fig. [TTT j-k)]. Note that the type of the elec- 
tron localization observed here in the second-excited state [Fig. [TT] (g-i)] corresponds to the 
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two-electron ground state for the larger value of d = 3.25. 

Although in all the low-energy states discussed here for d > 3.5 the electrons occupy 
different rings, in fact the type of the interring correlation between the electron positions 
varies significantly from state to state. Fig. [T2] shows the pair correlation function for spin 
singlet states calculated at d = 3.75 and B = 0.51 T for the ground state [Fig. [TWa-c)]. and 
the excited states [Fig. [T2](d-1)]. In each of the three columns we assumed that the electron 
in the right ring is localized in a different position - marked by the dot in Fig. [121 In the 
ground-state the electron in the right ring tends to stay far from the left ring and the angular 
reaction of the electron in the right ring on the actual position of the the electron in the left 
ring is weak [Fig. [T^ a-c)]. Only a slight rotation of the electron probability distribution in 
the left ring is observed in Fig. [T^ b) for the electron in the left ring localized at = 7r/2 
angle. For the ground state the overall PCF values within the right ring are reduced when 
the electron in the left ring is localized closer to the right ring [see Fig. [T2]^a-c) - contour 
scale is kept the same in all the plots]. In the third excited state an opposite tendency is 
observed [see Fig. [T2](j-1)]: the overall value of the electron probability distribution increases 
when the electron in the left ring approaches the right ring. In this state the minimal PCF 
value within the right ring is found always near = 0. In contrast to the ground state 
and the third excited state a strong correlation in the angular positions of the electrons is 
observed for the first [Fig. [TWd-f)] and second [Fig. [TWg-i)] excited states. 

In the first excited state the electron in the right ring stays away of the left ring when 
the electron in the left ring is at = [Fig. [T2l(f)]. which is natural for repulsing particles. 
However, when the electron in the left ring is localized above its center = 7r/2, the electron 
in the right ring is also most probable to be found above the center of its ring [Fig. [T2](f)] 
near = 7r/2, although one should rather expect a maximum of the probability density at 
the opposite side of the center. Also in the second excited state the correlation is somewhat 
different from what one might expect of the repulsing electrons. For instance when the 
electron in the left ring is localized at the closest distance to right ring the electron in the 
left ring tends to approach the left ring [Fig. [T2r g-i)]. 

In order to explain the observed features of the electron-electron correlation in the four- 
lowest energy singlet states one needs to consider the two-electron wave function. Below 
we present an approximate analysis for B ^ 0.5 T, when the total wave function is mainly 
constructed of the single-ring single-electron states of / = and / = 1 angular momenta 
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[see Fig. H], i.e. the single-electron ground state and the first excited state. Other energy 
levels lie much higher in the energy. The weak reaction of the electron in the right ring to 
the angular position of the electron in the left ring observed in the ground-state indicates 
that the two-electron wave function is nearly separable. The spatial part of the singlet 
ground-state wave function can be written in a following approximate non-normalized form 

^ = (i + Pi2)xKn)[/^(0i)-c/f(0i^ 

xx.(r2)[/^(02) + c/[(02)], (9) 

where is the electron exchange operator, '^-^ Xa is the radial wave function of the confine- 
ment within the ring a {a = I for the left and r for the right ring), and is the angular 
momentum m eigenstate for the electron localized in the ring a fm{<P) = exp{im(f)) / \/2ti. 
The argument of is the angular coordinate measured with respect to the center of ring a. 
In Eq. (j9]) the real valued c > introduces mixture of / = 1 to / = eigenstates. Mixture 
of these states is no longer rotationally invariant. For wave function ([9]) the electron in the 
left ring is localized preferentially near = vr and the electron in the right ring near = 0. 
For spatially separated xi Xr radial functions, the two-electron density is given by 



vl/|^ = (1 + Pi2)|xz(ri)xr(r2)r /o(0i) - c/l(0i 
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X|/o^(02)+c/[(02)|'. (10) 

For the first electron fixed within the left ring, Xri^i) ~^ and P12 may be skipped of the 
above formula since the permutated term vanishes.— We are left with the density which is 
a separable product of single-electron densities whose angular dependence on 02 becomes 
independent of 0i, although the overall pair correlation function value decreases when 0i 
approaches (the right ring), in agreement with Fig. [12]^a-c). 

Now let us turn our attention to the third excited state. In this state the electrons tend 
to occupy the area of the contact between the rings and the PCF shows a weak angular 
correlation between the electrons [Fig. [TWj-l)] like in the ground state. The wave function 
which produces these properties differs from Eq. (Q by the sign of c 

= a + Pi2)\Xi{ri)Xr{r2)\'\f!>{<Pi)+cf [{<!>,) 

X|/o(02)-c/[(02)r. (11) 

for which the electron localization angles are inverted (shifted by tt) with respect to the 
ground state. 
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In the ground state and in the third excited state the / = and / = 1 single-electron 
angular eigenstates are mixed within each of the rings and the resulting probability densities 
are (nearly) separable. In the first and second excited states the angular momenta / = and 
/ = 1 eigenstates contribute in a different manner. In these states the angular momentum 
of the electron in one of the rings is / = while the angular momentum of the electron in 
other ring is / = 1. This case is described by the following wave function 

^ = (1 + P12) [xK^i)/o(0i)Xr(r2)/[(02) ± Xi{ri)fi{<Pi)Xr{r2)m2)] (12) 

(for brevity we concentrate on the spin-singlets only). At B = 0.51 T for d = 3.75 the first 
(second) excited singlet is of the odd (even) parity, which corresponds to the '+' ('— ') sign 
in the above formula. In the absence of the overlap between the single-ring wave functions 
the two-electron density for the first electron fixed in the left ring is given by (up to the 
normalization constant) 

1^1' - IX/(ri)x..(r2)|' (1 ± cos(<^i - ^2)) , (13) 

where 0i is measured with respect to the center of the left ring and 02 with respect to the 
center of the right ring. Formula (fT3!) indicates that in the first excited state (of the odd 
parity) both the electrons are most probably localized at the same angle 0i = 02, and in 
the second excited state (of the even parity) the electrons are localized at opposite angles 
01 = 02 + TT, which explains the behavior observed in Fig. [121 g-i). At a close inspection in 
Fig- [TWh) one notices that the exact angular position 02 of the maximum in the right ring 
deviates off 0i a little bit, which is due to the contribution of the higher angular momenta 
to the wave functions. We see that in contrast to the ground-state and the third excited 
state the angular correlations between the electrons are strong in the first and second excited 
states. According to the presented arguments the symmetrization of the wave function by 
(1 + P12) operator [as well as the antisjTiimetrization by (1 — P12)] does not influence the 
correlated properties of the system when the electrons are spatially separated. Then, i) 
the singlet and triplet energy levels become degenerate and ii) the type of the correlation 
depends only on the form of the original wave function prior to symmetrization by (1 ±Pi2). 
The wave function of the first and second excited states that is symmetrized in [Eq. f[T^ ] 
is essentially entangled, hence the strong interring angular correlation. For the ground- 
state and the third excited state with negligible angular correlations the symmetrization 
was performed on a strictly separable product of single-electron wave functions. 
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FIG. 12: Pair correlation function plots for spin singlets with the position of the electron in the left 
ring fixed in the position marked by the dot. The distance parameter is d = 3.75 and the magnetic 
field B = 0.51 T. Lower rows of plots correspond to lower energies. 




FIG. 13: Three-electron spectra for a pair of identical rings. Solid and dashed curves show the 
even and odd parity energy levels. Black and red lines correspond to the total spin S = 1/2 and 
S = 3/2. 

D. Three electrons in a pair of identical rings 

The spectrum for three electrons at wide interring barrier is given in Fig. [T^ a). In 
the ground-state energy level at d = 3.75 we notice crossings of energy levels related to 
symmetry transformations which were also present in the single-electron spectrum but which 
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were absent for two-electrons. We notice that the "period" of the ground-state energy level 
oscillation is halved with respect to the single ring case [cf. Fig. [6]. In the low-energy 
spectrum we have two electrons in one ring and a single electron in the other. Halving of 
the oscillation period results of the fractional Aharonov-Bohm effect^ occurring for few- 
electron states confined in a single ring. The ground-state oscillations are due to the spin- 
transformations of the electron pair in a single ring. The angular momentum of the two- 
electron subsystem is not a good quantum number due to the perturbation of the two- 
electron ring by the potential of the electron in the right ring. For 5 < 0.5 T the ground 
state is two-fold degenerate: with respect to the parity and corresponds to S* = 1/2. For 
low magnetic fields the spins of the electrons in the single ring are opposite and compensate, 
so the total spin can be identified with the spin of the solitary electron in the other ring. 
For B > 0.5 T the ground state of the two-electron subsystem is the spin-triplet. Since the 
spin of the solitary electron may have an arbitrary orientation the ground-state becomes 
four-fold degenerate: with respect to both parity and the both allowed total spin quantum 
number values S = 1/2 and 3/2. For smaller d the degeneracies are lifted [see Fig. [T3l(b-d)]. 
Near the ground-state we notice a characteristic oscillation of the ground-state symmetry 
[Fig. [T3](c.d)]. The ground-state is mostly either the odd-parity low-spin S* = 1/2 state or 
the odd-parity spin-polarized state S = 3/2. When the energy order of these two energy 
levels changes they become nearly degenerate with a third state: the even-parity low-spin 
state. The even-parity spin-polarized {S = 3/2) state runs much higher in the energy. This 
sequence of the ground-state spin- and parity-symmetry oscillations is also characteristic 
to three-electron circular— i^i^ dots [for circular dots states of angular momentum quantum 
number which is even (odd) integer are of the even (odd) parity] as well as to elliptic^ 
quantum dots and double quantum dots^^ containing three electrons. For elliptic quantum 
dots the near degeneracy of the three energy levels that we obtain here for discrete values of 
the magnetic field indicates that the deformation of the confinement potential with respect 
to the circular limit is strong (see the discussion given in Ref.— for elliptic dots). For less 
strong deformation the even parity low-spin state becomes a ground-state for some narrow 
but distinct magnetic field ranges.—"^ 
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FIG. 14: Ground-state charge density for three electrons in identical rings (d-f) and for the right 
ring deeper by 0.1 meV (a-c) at the interring spacing of d = 3.75. 

E. Oscillations of the three-electron charge density with the magnetic field 

The results for the charge density of one- and two- electron systems presented above 
exhibited a distinct magnetic field dependence. The results for three electrons at d = 3.75 
presented in Fig. [MT d-f) indicate that the dependence on the magnetic field is significantly 
weaker. For ideally symmetric pair of rings the charge density of both even and odd parity 
three-electron eigenstates is distributed equally between the rings with 1.5 electron charge 
per ring on average. A classical distribution with one electron in one ring and two in the 
other can only occur when the symmetry is lifted. In the absence of the tunnel coupling a 
potential well difference of 0.1 meV is enough to obtain the integer distribution of electrons 
between the rings. The charge density for = 3.75 and the right ring deeper by 0.1 meV is 
shown in Fig. [TWa-c). The dependence of the single-electron charge localized in the left ring 
on the magnetic field is clear, and the reaction of the two-electron density in the right ring 
is weaker. Note that in this plot the deviation of the single-electron density off the circular 
symmetry occurs for B = 0.91 T, for which the two-electron density is closest to circular, 
which leads to a compensation of the magnetic oscillation of the density for identical rings. 

In order to quantify deviations of the single-ring electron density off the circular symmetry 
we calculated the parameter 

D =< {x - X)^ > - < {y -Yf >, (14) 

where X and Y are coordinates of the ring and the average is calculated over the charge 
density of separate rings. The result is displayed in Fig. [16] with the blue curve for the 
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FIG. 15: Pair correlation function for the ground-state at = and d = 3.75 for the right ring 
deeper by 0.1 meV than the left one. Left ring contains a single electron and the right one two 
electrons. A position of one of the electrons in the right ring is fixed and marked by the dot. 
Contours of the ring area are also shown. 




0.4 0.8 1.2 1.6 

Bm 



0.4 0.8 1.2 1.6 2 
B[T| 



FIG. 16: (a) Deviations o the single-ring electron densities off the circular symmetry [Eq. (114p ] in 
function of the magnetic field for the system of three electrons in double ring structure with the 
right ring deeper by 0.1 meV for d = 3.75. Blue curve corresponds to the deeper ring occupied 
by two electrons and the black curve to the shallower ring which contain a single electron, (b) 
Derivative of D over B. a_B is the Bohr radius. 

two-electron density of the deeper ring and with the black curve for the single-electron 
density of the shallower ring. The single-electron parameter is larger, varies more strongly 
with the magnetic field and its variation is continuous in contrast to the two-electron value, 
which has cusps when the ground-state spin state of the two-electron subsystem changes. 
Moreover, the single-electron parameter is positive which indicates that the deformation 
occurs rather in the horizontal (x) direction, while the two-electron deviation occurs mainly 
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in the vertical direction (y). The absolute value of the deviation of the two-electron density 
from circular is the largest at the symmetry transformations. As noted in context of Fig. 
[m the largest deformation of the single-electron density corresponds to a weak deformation 
of the two-electron density. 

When the two-electron subsystem changes its spin state the ground state charge density 
in the deeper ring is modified in a discontinuous manner. Due to the Coulomb coupling 
this change may influence the electron density of the shallower ring, which might react 
in a discontinuous manner. In order to quantify this reaction we calculated derivatives of 
the D parameter with respect to the magnetic field. The derivative of the two electron 
parameter is discontinuous at the cusps [see Fig. [T^ b)]. We notice that the single-electron 
parameter exhibits an irregular structure when the two-electron density changes. However, 
this structure is not very pronounced. This result along with the presented above reference 
calculations indicate that the charge density of each ring is quite indifferent to the actual 
form of the charge distribution in the other ring. Additional argument for that conclusion 
is provided by the pair correlation function. 

The pair correlation function plot for three electron system with two electrons in the right 
ring is shown in Fig. [15] for the ground-state at zero magnetic field. We fix a position of 
one of the electrons in the right ring and mark it with the dot in the figure. We notice that 
the angular distribution of the solitary electron in the left ring is unaffected by the position 
of the fixed electron. Only the electron in the same (right) ring reacts to the variation 
of the fixed electron position [compare Fig. [T5l(a) and fT5l(b)]. Note, that for the electron 
fixed below the center of the right ring [Fig. [T6](a)] the other electron in the right ring is 
not exactly on the other side of the center. It tends to avoid the left part of the ring and 
the position of the maximum is localized below the angle = 7i/2. The result of Fig. [15] 
indicates that in rings separated by a barrier that is thick enough to prevent the interring 
tunneling the electrons are mainly coupled as entire charge densities and do not react on 
their actual position. A strong angular correlation is only observed for the electron within 
the same ring. 
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F. Magnetization 



Fig. [T7] shows the magnetization M = — |^ of the pair of identical rings filled by 1 , 2 and 
3 electrons with different spacing parameters d. For a single circular ring the magnetization 
is discontinuous at the angular momentum transitions. The interring tunnel coupling makes 
the single-electron magnetization smooth and continuous [see black curves in Fig. [T71(a-d)]. 
The oscillation of M with B are extinguished for the coalesced rings, when the Aharonov- 
Bohm oscillations occurs only in excited states and not in the ground-state which is localized 
in the quantum dot formed at the contact of the rings. 

For two electrons the magnetization is a smooth function of the magnetic field unless 
the rings are close enough to form a single-structure. Two-electron ground-state symmetry 
transformations which produce the cusps in the magnetization occur only for the coalesced 
rings {d = 3). For three electrons the symmetry transformations are present for any d. For 
large interring distance they correspond to the spin transformations of the subsystem of 
two electrons confined within the same ring. The spin transformations of the three electron 
system occur also when the rings form a single coalesced structure. 




FIG. 17: Magnetization of a pair of identical rings filled with 1, 2 or 3 electrons plotted with black, 
red and blue curves, respectively. 



G. Charging properties 

For a structure embedded in a charge tunable device the charging of the double ring by 
N — th electron occurs when the chemical potential = En — E^^i (where En is the 
ground-state energy of confined electrons) is aligned with the Fermi level of the electron 
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reservoir. In capacitance spectroscopy^ the cusps of the charging hnes in function of the 
magnetic field indicate symmetry transformations of the confined system. 

Fig. [TH] shows the calculated chemical potentials in function of the magnetic field. Chem- 
ical potentials of 1 and 2 electrons were shifted up on the energy scale for clarity (the amount 
of the shift is given in the figure). In the absence of the tunnel coupling [Fig. [TBlfa.b)] the 
chemical potential of two electrons is a smooth function of the magnetic field with exception 
to the V shaped cusps that result of the single-electron angular momentum transitions that 
produce the A shaped cusps on the single-electron chemical potential. Since the two-electron 
ground-state energy is smooth as a function of B all the cusps of the three electron system 
have the A shape and result of the symmetry transformations of the three electron system. 
The parabolic minimum of /is between the A cusps is due to the smooth maximum of E2 
that occurs at the avoided crossings near odd multiples of flux quanta (see Fig. [9]). Due 
to the fractional Aharonov-Bohm oscillation of the two electron subsystem we find two A 
cusps in /is for a single one in /xi. 

For d = 3.25 the A cusps of the single-electron chemical potential are smoothed out and 
in consequence the V cusps of /i2 disappear [Fig. [TSlfc)]. The ground-state avoided crossing 
for the two-electrons becomes narrower, so fi2 acquires a distinct maximum near 1 T. A A 
cusp in /i2 and consequently a V cusp in /is occur for the coalesced rings [Fig. [T8l(d)]. The 
presented results indicate that the strength of the interring coupling may be deduced of the 
charging experiments. 




FIG. 18: Chemical potential of 1, 2 and 3 electrons confined in the double ring plotted with black, 
red and blue curves, respectively. 
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H. Separability of the system at large interring barrier 

In the absence of the tunnel couphng the rings are coupled only electrostatically. In order 
to answer the question to which extent the distinct rings may be treated as separable we 
performed a reference calculations in which we considered a single ring perturbed by the 
Coulomb potential of a classical point charge localized in the other ring. In the reference 
calculation for the two-electron system we assumed that in the left (right) ring a classical 
point-charge electron is localized at = vr (0 = 0) and calculated the single-electron spec- 
trum for the confinement potential of the right (left) ring. The two-electron spectrum of a 
ring couple was then estimated by the sum of single-electron spectra of the left and right 
rings. 

Comparison of the obtained result with the exact two-electron spectrum is given in Fig. 
[T9l The reference calculation was shifted down by 0.41 meV to coincide with the exact two- 
electron energy for B = 0. We see that the magnetic- field dependence of the ground-state 
energy and the nearly degenerate first excited energy level is quite accurately described by 
the ansatz model. In fact the exact first excited energy level is not exactly degenerate and 
there are two energy levels for states of different parities. The reference calculation in which 
the rings are treated separately overlooks this splitting. Although the width of the exact 
avoided crossing near i? = 1 T is quite accurately described by the reference calculation 
the increase of the exact ground-state energy in the center of the avoided crossing is smaller 
than in the reference calculation. The higher part of the exact spectrum deviates off the 
reference calculation significantly. 

The reference spectrum for three electrons was calculated as a sum of two spectra 1) 
the single-electron confined within the left ring which is perturbed by two classical charges 
localized in the right ring at positions marked by the crosses in the blue inset to Fig. [20] and 
2) the two-electron spectrum of the right ring with the classical point charge localized in the 
left ring in the position of the cross in the black inset to Fig. [201 Blue curve in Fig. [20] shows 
the single-electron ground-state, black and red curves correspond to the lowest-energy singlet 
and triplet states of the two-electron ring, respectively. For a single circular ring the confined 
electron system undergoes momentum transitions as function of the magnetic field. For the 
two-electrons the transitions are accompanied by the spin transformations. When the two 
rings - one containing a single electron and the other two electrons - become electrostatically 
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FIG. 19: The exact two-electron spectrum (lines) and the reference spectrum (dots) obtained as 
a sum of single-electron spectra for a single ring perturbed by a classical charge within the other 
ring. The reference spectrum was shifted down by 0.41 meV on the energy scale. The crosses in 
the inset show the assumed positions of the classical charges. 




FIG. 20: Blue line shows the energy of the single electron (shifted down by -44 meV) within the 
left ring with the Coulomb potential of two classical point charges localized as indicated in the 
blue inset to the figure. Black and red curves show the lowest-energy singlet and triplet for two 
electrons in the right ring with a classical point charge localized in the left ring in the position 
marked in the black inset. Distance parameter between the rings is d = 3.75 and the depth of the 
right ring is increased by 0.1 meV. 

coupled the angular momentum transformations disappear. However, the spin ground-state 
transitions are still present in the two-electron subsystem. 

The reference spectrum for the three-electron system calculated using the two separate 
calculations is compared to the exact spectrum in Fig. [211 In the ground-state both calcu- 
lations agree quite accurately, but differences are observed in the excited states. 
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FIG. 21: The exact three-electron spectrum (hnes) and the reference calculation (dots). The 
reference spectrum is a sum of a single-electron spectrum in the left ring perturbed by two classical 
charges and the two-electron spectrum of the right ring with a classical charge localized in the left 
ring (see the inset to this figure and to Fig. I2U|) . Black color corresponds to the low-spin states 
and red color to the high spin states. The reference spectrum was shifted down by 1.41 meV. 

I. Electron-electron correlation in the coalesced double rings 

The pair correlation plots for the three electron system presented so far were limited to 
large interring barrier. Let us look at the correlation when the rings form a single structure. 
The pair correlation function plots for d = 3 and B = are given in Fig. [22] for the low-spin 
states. Near the ground state one of the electrons occupy the contact area between the rings 
and the two others are localized at the left and right ends of the double structure. In the 
plots we present results for two different locations of the electron in the left ring (marked by 
crosses in Fig [221) . The pair correlation for the lowest-energy odd parity state Fig. [22]^a,b) 
and the first excited state Fig. [22]^c,d) are similar. Both these states have similar energies 
(see the spectrum in Fig. [T3l(d)]. Only the electron localized at the contact between the 
rings reacts to the position of the electron in the left ring, and this reaction is not very 
pronounced [see Fig. [22](a-d)]. A stronger reaction is observed in the three next excited 
states [Fig. [22l(e-j)]. which correspond to distinctly higher energies [Fig. [T3l(d)]. In the 
highest energy state of the presented set the "quantum dot" formed at the contact is empty 
and the correlation properties [Fig. [22](i.j)] are similar to the ground-state at large interring 
barrier [cf. Fig. [15]. 

Results presented above for wide interring barrier indicated that in the ground-state the 
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actual electron positions are correlated only within the same ring. Fig. [22] demonstrates that 
for rings forming a single structure the electron-electron correlation in the ground state still 
has only a short-range character. In the studied case each of the three electrons occupy a 
different spatial location and is quite indifferent to the actual positions of the other electrons 
within their charge islands. 
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FIG. 22: Pair correlation function for three electrons in coalesced double ring structure d = 3 
at = 0. Plots (a,b) correspond to the ground-state, and higher rows to subsequent low spin 
(S = 1/2) energy levels. The two columns correspond to two locations of the fixed position 
electron which is marked by the cross in the figure. 



J. Rings of different radii 



Results presented so far were obtained under assumption that both the rings have the 
same size. Let us now consider two rings that possess slightly different radii. Namely, we 
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assume that the radius of the right ring is increased by 5 percent. The calculated single- 
electron spectra are displayed in Fig. [231 

The energy of the lowest zero angular momentum state in a single ring increases with 
its increasing radius (see Fig. [2^ . For that reason in the absence of the magnetic field the 
zero-angular momentum ground-state is localized in the smaller ring [energy levels marked 
by S in Fig. [23]^a,b)] and the first excited state is localized in the larger ring (energy levels 
marked by L in Fig. [23]) . On the contrary, due to the centrifugal potential at zero magnetic 
field the states of the non-zero angular momentum have lower energy when localized in 
the larger ring [see Fig. [231(a)]. When the magnetic field is switched on the energy levels 
that originate of the positive angular momenta single-ring states go down in the energy 
due to the interaction of their paramagnetic dipole moments with the external field [see the 
last term in the Hamiltonian ([HI)] and they replace the zero-angular momentum low-field 
ground-state when B is high enough. Consequently, for B exceeding 1 T the ground-state 
becomes localized in the larger ring. For 5 > 1.2 T the ground state (in the larger ring) 
and the first excited state (localized in the smaller ring) start to increase in the energy. This 
increase results of the diamagnetic shift [related to the second term at right of (Eq. [HI)]. It 
is a general rule that the diamagnetic shift is stronger for more delocalized states (covering 
areas of larger p). In consequence one obtains another change of the electron localization: 
above B = 1.8 T the ground state is localized again in the smaller ring. The ground- 
state transition to the larger ring occurs again when the states originating of higher angular 
momenta become ground-states. For d = 3.75 the change in the ground-state localization 
occurs through crossing of energy levels [see Fig. [231(a)]. Already for d = 3.5 tunnel-related 
avoided crossings are obtained [see Fig. [23r b)]. For the strong coupling case of d = 3.25 the 
spectrum resembles the one of the identical rings [cf. Fig. [6]^c)]. 

The shifts of the electron density with the magnetic field between the rings are illustrated 
in Fig. [25l for d = 3.5. For low magnetic field the electron in both the ground-state and 
the first-excited state is present in both the rings [Fig. [25t^a) for B = 0.5]. At B = 1 
T we observe an avoided crossing which results of the change of the order of the energy 
levels corresponding to different single-ring angular momenta. For B < 1 T {B > 1 T) the 
electron in the ground-state occupies preferentially the smaller (larger) ring [see Fig. [25]^a) 
and [25r b)]. As the magnetic field increases above 1.25 T the ground-state presence of the 
electron in the left (smaller) ring is enhanced. Equal electron distribution is obtained near 
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the center of the avoided crossing for B ~ 1.75 T [Fig. [251(c)] and for B = 2.51 T the 
ground-state is totally localized in the smaller ring. Near B = 2.8 T in [Fig. [23r b)] we see a 
change of the order of the energy levels in the ground-state. This is a very narrow avoided 
crossing (which in the figure scale appears as a level crossing). After this avoided crossing 
the electron in the ground-state is localized in the larger ring again. Equal distribution of 
the ground-state electron between the two rings is found near 3.61 T. The change of the 
order of the two energy levels in the ground-state that occurs here is also accompanied by 
a narrow avoided crossing. For the magnetic field increased additionally by only 0.2 T the 
interring tunnel coupling is broken and the ground-state and the first-excited state occupy 
a single ring. We see that at higher magnetic fields - which favor a stronger localization 
of the confined states - the avoided crossing related to electron transfer between the rings 
become narrower both on the magnetic field and energy scales (compare Fig. [25r a-c) and 

Fig. ESKd-f)]. 
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FIG. 23: Single-electron spectrum for a pair of rings of identical depth but with the radius of the 
right ring increased by 5%. In (a) and (b) states localized mostly within the smaller (larger) ring 
are labeled by S (L). In (c) electron in the ground state is nearly equally distributed between the 
rings. 

To summarize the shifts of the charge in the single-electron problem of a single-ring in 
Fig. [26] we plotted the ground-state charge localized in the left ring for various distance 
parameters. For d = 3.75 the charge in the left ring has a step-like dependence on B and 
is either very close to 1 or very close to zero. The oscillation of the charge for d = 3.5 
becomes less abrupt. For the strong coupling case d = 3.25 the charge localized in the left 
ring oscillates around 0.5e. In all the cases the oscillations of the charge localized in the left 
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FIG. 24: Ground-state energy in function of the ring radius for a single electron within a single 
ring. 




(c) 





(f) 

FIG. 25: Cross section of the confinement potential (grey lines) and the charge density of the single 
confined electron in the lowest-energy state (blue curve) and the first excited state (red curve) for 
a pair of quantum rings with the radius of the right ring larger by 5% and the interring distance 
parameter d = 3.5. Panels correspond to different values of the magnetic field given in the figure. 

ring become more pronounced as the magnetic field grows, which is related to the reduction 
of the tunnel coupling with B. 



31 



FIG. 26: Charge localized in the smaller ring for ground state of the single-electron (a) two-electrons 
(b) and three-electrons (c) for the ring couple with the radius of the right ring larger by 5%. Red, 
black and blue lines correspond to d = 3.75, d = 3.5 and d = 3.25, respectively. 

For two-electrons the Coulomb repulsion segregates the electrons between the different 
rings. For d = 3.75 and d = 3.5 the charge confined in each of the rings is equal to le [see 
Fig. [26r b)]. Only for d = 3.25 deviations off the equal distribution of the rings occur. We 
find that these deviations are accompanied by lifting the singlet-triplet degeneracies in the 
spectrum. The magnetic field leads eventually to attenuation of the tunnel coupling and 
segregation of the carriers between the two dots. 

In the system of three electrons [see Fig. [26]^c)] for large interring barrier {d = 3.75) 
two-electrons occupy the larger ring leaving a single electron in the smaller one. For d = 3.5 
the charge accumulated in the smaller ring exceeds one elementary charge for low magnetic 
fields. For d = 3.25 the electron distribution between the rings becomes closer to 3e/2 
per each ring. Oscillations of the charge in the smaller ring for d = 3.5 and d = 3.25 are 
associated with the avoided crossings that appear in the low-spin part of the spectrum as 
well as with the ground-state spin transitions. The latter produce cusps on the plotted 
curves [see Fig. 26(c)]. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the system of one, two and three electrons in a planar double ring 
structure considering both the tunnel and the electrostatic coupling between the rings using 
a Gaussian functions mesh technique. 
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The presented results indicate that in the system of few electrons there is a distinct 
competition of tunnel coupling with 1) the magnetic field which enhances localization of the 
occupied states within a single ring 2) the electron-electron interaction which strengthens the 
interring barrier and 3) asymmetry effects of the double ring structure which favor electron 
localization in one of the rings. 

We find that the system undergoes symmetry transitions of the parity and spin in function 
of the magnetic field, and that these transitions vary strongly with both N and d. At large 
interring barrier (d) ior N — 1 we find angular momentum transitions for states localized in 
a single ring. These transitions are absent in the system of two electrons since their mutual 
interaction perturbs the rotational symmetry of each of the rings. N = 2 ground-state at 
large d is degenerate with respect to the spin due to the perfect separation of the electrons. 
For = 3 at large d the symmetry transformations are present, and result of the spin 
transitions within the two-electron subsystem confined in one of the rings. When the rings 
are identical at large d the ground state is degenerate with respect to the parity for any N. 

In the opposite limit of strong interring coupling and coalesced rings we find lifting of the 
even-odd degeneracy for all the studied electron numbers. The symmetry transformations 
in function of the magnetic field vanish for the single electron, which occupies the binding 
orbital and in the limit of coalesced rings - the quantum dot formed at the contact of the 
rings. On the contrary, for two electrons the symmetry transformations occur only when 
the rings form a single structure. For three electrons the symmetry transformations are 
present at any d. At large d they are only related to the spin transitions of the two-electron 
subsystem and in the strong coupling limit they involve both parity and spin. For d = 3 the 
sequence of the parity and spin ground-state symmetries in function of the magnetic field is 
identical with the one found for strongly deformed elliptic quantum dots. 

We demonstrated that due to the strong dependence of the symmetry transformations on 
both A^ and d, the confined charge as well as the distance between the rings should be readily 
accessible in both single-electron charging and magnetization measurements performed in 
function of the magnetic field. In particular the strength of the coupling has an opposite 
impact on the chemical potentials of one and two confined electrons. For week tunnel 
coupling /ii has cusps in function of the field while H2 is a smooth function of B. In the 
strong coupling limit the dependence is inverted: chemical potential of the single electron 
becomes smooth and the cusps appear on ii2- 
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The asymmetry in the depth of the confinement potential favors locahzation of electrons 
in the deeper potential well. Less obvious is the effect of the asymmetry in the size of 
the rings. For rings of different radii the magnetic field leads to switching of the ground- 
state localization from the larger to the smaller rings. The oscillations of the ground-state 
localization become more abrupt for stronger magnetic fields due to the reduction of the 
tunnel coupling. Both types of the asymmetry - in the depth and in the size of the ring 
- for large values oi d oi B - break the tunnel coupling between the rings, by tending to 
localize the electron in one of the rings - favoring the atomic type of localization over the 
molecular (extended) one. The asymmetry in the double ring potential lifts the degeneracy 
of the even and odd parity energy levels that is observed at large d for any A^. We find 
that when d is decreased, the asymmetry effects are reduced. In particular the electron 
charge between the rings becomes more evenly distributed for the strongly tunnel-coupled 
asymmetric structures. 



V. APPENDIX 



In this appendix we explain how the matrix elements for the Coulomb interaction 

Vabcd = (0a(ri)0b(r2)| — |0c(ri)0d(r2)), (15) 

12 

are integrated. In the above expression the single electron wavefunctions are replaced by 
their linear combinations ([3]): 

Va,ai= E <6*CfcC?,(/,(ri)/,(r2)| — |/,(ri)/Kr2)) (16) 

i,j,k,l=l ^12 

An interaction integral which appears in (fT6|l : 

Q,ki = (/.(ri)/,(r2)| — |/fc(rO/Kr2)) = (17) 

12 



/ dh, [ rfV/;(ri)/;(r2)— /,(rO/Kr2) 



ri2 

can be calculated in following way. First, we should substitute if place of /j, /j, /fc, // and 
their invert Fourier transform. Next step is to integrate over ri and Y2 variables. After 
these we get the formula on Cijki- 

Cijki = B j rf^ki exp - ^) exp ■ r) exp (^^k ■ r) (18) 
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where 



R 



Rjfc 

Rj7 

7 



R{cos P, sin P) 
iVji - Vik, Xik -Xji) 
r(cosa, sin a) 

Ri ~ Rfc 

Rj — R« 
Rj + Rfc 

2 

Rj + R/ 
2 

2/r 



Ob is Bohr radius, 
and 



B 



na 



exp 



X exp (^i^{-YikXik - K,7a;j7 + XikVik + Xjiyji)^ 



X exp 



rfk + r]i 
4a2 



The expression [TSl can be further simphfied by doing transformation the integral form 
sian to the cyhndrical coordinates and then by integrating over the k variable: 



ijkl 



B 



2a Jo 

The function A{(p) is defined as below: 



dip exp 



2a2 



erf 



\V2aJ 



A{ip) = —iRcos{(p — P) 
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r cos(v9 — a) 



where 7 = — |f • The value of integral ( 1271) is calculated numerically. 
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